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When searching for deviations of statistical isotropy in CMB, a popular strategy is to write the 
two-point correlation function (2pcf) as the most general function of four spherical angles (i.e., two 
unit vectors) in the celestial sphere. Then, using a basis of bipolar spherical harmonics, statistical 
anisotropy will show up if and only if any coefficient of the expansion with non-trivial bipolar 
momentum is detected - although this detection will not in general elucidate the origin of the 
anisotropy. In this work we show that two new sets of four angles and basis functions exist which 
completely specifies the 2pcf, while, at the same time, offering a possible geometrical interpretation 
of the mechanisms generating the signal. Since the coefficients of these expansions are zero if and 
only if isotropy holds, they act as a simple and geometrically motivated null test of statistical 
isotropy, with the advantage of allowing cosmic variance to be controlled in a systematic way. We 
report the results of the application of these null tests to the latest temperature data released by 
the Planck collaboration. 


I. INTRODUCTION 

Given the unprecedented limits on cosmological pa¬ 
rameters achieved with Planck data releases mm , an 
important follow-up question is whether the same data 
contain traces of physics beyond the standard ACDM 
model. While deep field surveys Ena aim to unveil the 
specific nature of dark matter and dark energy, and thus 
of the energy content of our universe, CMB observations 
are special in the sense that they provide a unique win¬ 
dow to both the physics of the early inflationary universe 
and of its global shape, i.e., its geometry and topology E3- 

m- 

From the statistical point of view, finding evidences of 
new features of inflationary physics or of the shape of the 
universe usually translates into cosmological detections of 
non-Gaussianity and statistical anisotropy, respectively 
(but not necessarily [TO]). As it turns out, however, ob¬ 
servational bounds from CMB on non-Gaussianity |11) - 
Ea did not allow us to discriminate inflationary mod¬ 
els, and the hope now lies on the possibility that future 
measurements of polarization R-modes of CMB m will 
better elucidate the physics of the early universe. The 
physics describing the global topology and geometry of 
the universe, on the other hand, is not only much less 
constrained by CMB [5] , but is also equally fundamental 
to the beyond-ACDM program. In this regard, the exis¬ 
tence of statistical anomalies at the largest CMB angles 
(see Ref. [15] for a comprehensive review) can be opti- 
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mistically seen as an indication of spatial anisotropy at 
the horizon scales m, although the conservative minded 
would also remind us of possible unaccounted systematic 
effects m, or even the less exciting case of statistical 
flukes p~8] . 

Regardless of the final words on CMB anomalies, how¬ 
ever, one is rightfully justified to question the validity of 
the statistically isotropic scenario, given its deep connec¬ 
tions with the symmetry hypothesis about our universe. 
Thus motivated, this paper addresses the question of how 
to constrain deviations from statistical isotropy in a ge¬ 
ometrically meaningful and model-independent way. 

The general recipe for describing the statistics of a 
Gaussian and statistically isotropic CMB map is straight¬ 
forward. Assuming that the geometry of the universe is 
everywhere rotationally invariant, all we need to do is to 
compute the correlation between the temperature fluctu¬ 
ations AT at directions n and n ', given by 

C{h, n') = (AT(n)AT(h '))\ cos = C\0 ). (1) 

However, if there are departures from statistical isotropy, 
either of systematic HD, astrophysical HMD or cosmo¬ 
logical H3I22H23] origin, the two-point correlation func¬ 
tion (2pcf) will depend on other angles relating n and 
n'; in this case, C(6) no longer exhausts the statistics of 
the universe. Thus, if we want to go beyond the ACDM 
framework, the central question is how to parameterize 
deviations from Eq. ([l]) in a meaningful, and hopefully 
practical, way. 

The angular correlation function has come under con¬ 
siderable scrutiny. It was first noticed by the COBE- 
DMR team [25] that C{0) was unexpectedly close to zero 
for 0 > 60°. This lack of correlations was confirmed 
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by the Wilkinson Microwave Anisotropy Probe (WMAP) 
team in their analysis of their first year of data [125] , 
Though WMAP claimed to have greatly reduced signif¬ 
icance in future releases, Copi et al. showed m that 
in fact this absence of large-angle correlations persisted 
on the sky outside the galaxy in the third year release, 
and in all subsequent releases, including the first-year 
Planck release pzSHBU] . Those findings have since been 
confirmed by others EH 1323, although no satisfactory ex¬ 
planation exists. It has been suggested |33ft35] that one 
might be able to test whether the anomalous vanishing of 
the temperature-temperature correlation function is due 
to new physics or just a statistical anomaly by examining 
other two-point correlation functions (eg. the EMnode- 
A-mode correlation function) on similarly large angles. 

A common strategy for parameterizing deviations of 
isotropy in the CMB is to use a complete set of ba¬ 
sis functions to perform a multipolar expansion of the 
2pcf. Since the latter is defined by the product of two 
functions on the CMB sphere, we could simply use the 
product of two independent spherical harmonics as such 
basis. Instead, it has become a standard practice to 
use a basis of total angular momentum eigenfunctions, 
also known as bipolar spherical harmonics ESI) to do 
the expansion. Besides sharing most of the mathemat¬ 
ical properties of the standard spherical harmonics, the 
advantage of the bipolar harmonics is that they encode 
deviations of isotropy in the total angular momentum of 
the coefficients of the expansion. Thus, any measure¬ 
ment of a multipolar coefficient with a non-trivial total 
angular momentum eigenvalue is an indication of statis¬ 
tical anisotropy. This program was introduced by Hajian 
and Souradeep [I37H39] , and has been fruitfully applied to 
CMB since then |201 PTOl PTlTj . 

Since the bipolar spherical harmonics form a basis for 
square-integrable functions on the Hilbert space where 
the 2pcf is defined, they offer a very general framework 
for studying the statistics of the CMB. However, it has 
some limitations, too. First, the multipolar coefficients of 
the bipolar expansion, while serving as a null indicator 
of anisotropy, do not provide a physical interpretation 
of the underlying signal straightforwardly, and for that 
one has to resort to other tools m- Second, in a more 
symmetric situation - whether real or expected - it is not 
clear how the degrees of freedom of the bipolar spherical 
harmonics could be combined to reduce cosmic variance. 

This paper is based on the idea that, given two unit 
vectors rooted at the origin of the CMB sphere, two new 
and unique geometrical objects can be formed: a great 
circle and (the boundary circle of) a spherical cap - or, if 
we include the interior of the sphere, a disc and a cone. 
Using the set of angles defined by these objects, one can 
introduce a complete set of basis functions that charac¬ 
terize the 2pcf, and whose multipolar coefficients offer a 
direct and geometrically motivated null test of statisti¬ 
cal isotropy. Moreover, since the angles used to represent 
the 2pcf have a clear geometrical interpretation, it may at 
times be physically well-motivated to integrate, average 


or marginalize over one of them. 

This work is organized as follows. After reviewing the 
basics of the bipolar spherical harmonic formalism in Sec. 
|II A| we introduce the new geometric representations of 
the 2pcf in Secs. |IIB| and |IIC[ where we also show how 
they recover the usual 2pcf in the isotropic limit. In 
Sec. ED we show how these new functions can be used 
to construct null geometrical tests of statistical isotropy. 
Finally, we present the results of null tests applied to the 
latest temperature maps released by the Planck team in 
Sec. |III A| We conclude in Sec. |IV| with some perspec¬ 
tives of future developments. 

II. GEOMETRIC REPRESENTATIONS 

We start by recalling the basic definitions and nota¬ 
tions used in CMB statistics. The temperature fluctua¬ 
tion field is a real function on the sphere, and can thus 
be expanded as 

A T(fi) = ^2 aimYimi.fi) , 

l,m 

where n = (0,<f>). In the canonical ACDM cosmological 
model, the eq m are realizations of a Gaussian-random 
and statistically independent variables. In this case the 
expectation value of the two-point correlation function is 

C (ni,n 2 ) =5Z<a; 1 m 1 < 2m2 )y) lTOl (ni)y ; * m2 (n 2 ). (2) 

h ,171l 

h,rri2 

The set of coefficients (ai imi a* 2m2 ) form the covariance 
matrix, and in the canonical statistically isotropic case 

( a Zimi a i 2m2 ) = ^-'h^hh^m 1 m2 ■ (3) 

Any non-zero off-diagonal term in the covariance matrix 
is a measure of statistical anisotropy. 1 
The 2pcf is symmetric by definition 

C(ni,n 2 ) = C[n 2 , ni). (4) 

As we shall see, this symmetry imposes restrictions on 
the eigenvalues of the eigenfunctions (and consequently 
on the multipolar coefficients) that we introduce below. 

A. Bipolar representation 

A convenient basis for expanding the 2pcf is given by 
the bipolar spherical harmonics (3BJj which are defined 


1 Conversely, a diagonal matrix does not imply isotropy, since Ci 1 
could depend on m\. Statistical isotropy is thus a strong con¬ 
dition requiring both independence among multipoles and the 
invariance of Ci 1 by rotations. 
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as the tensor product of two spherical harmonics 


®(n i,n 2 )= Y h (n i)®Y h (n 2 ) 


J LM 


where Yq is a shorthand notation for 
(Yq _q;...; Yq + q). In terms of this basis the 2pcf reads 

C(h 1 ,n 2 )= ( 5 ) 

L,M li ,l 2 

The coefficients are called the BipoSH spec¬ 

trum |37ll3)!)| . They are given by a quadratic combination 
of the a lm \ 

A hh = J2 C h™lhm 2 ( a lim i a l2 r n2 ) , (6) 

m 1 ,771 2 

where C, L ^ , are the Clebsch-Gordan coefficients. In 

1 3-222 1 * 2 ?™ 2 

the case pi) of statistical isotropy they reduce to 

A hh = 1 Cq 5 

so that any statistically significant detection of a non-zero 
with L > 0 is a sign of statistical anisotropy. This 
property makes the BipoSH coefficients a very convenient 
null test for statistical anisotropy. 


B. Anisotropies through conic modulations 

Given two unit vectors, hi = (xii<Ai) and f >-2 — 
(X 2 ; 4 * 2)1 rooted at the origin, they define a cone on the 
unit sphere, obtained by rotating those vectors about an 
axis collinear to hi 2 oc hi T- fi 2 . Each cone can be com¬ 
pletely described by three angles: the opening angle x of 
the cone, and two angles 0 and $ giving the orientation 
of the axis h 12 . The ordered pair of directions, ( 711 , 77 . 2 ), 
is fixed by a fourth angle 77 specifying their position on 
the circle bounding the intersection of the cone with the 
unit sphere 2 - see Fig. [l] Thus, instead of representing 
the degrees of freedom of the 2pcf by the usual spherical 
angles, as in Eq. ([ 2 ]), we can use the four angles defined 
by the cone: 


C(hi, h 2 ) = C(x, r], 0, <&). (7) 

These angles range over 

0 < 0 < 7T , 0<$<27T, 0<X< 7r ) 0 < 77 < 27T . 

Note that by construction x is the usual angle between 
hi and 77.2 

cos x = cos Xi cos X 2 + sin xi sinX 2 cos(<^i - 4> 2 ) ■ (8) 


2 We adopt the convention that rj is the angle between the arc of 
the great circle from hi to fi 2 and the arc of the great circle from 
—z to z through rq 2 . 



Figure 1: Two unit vectors in the CMB sphere will always 
define cones. The number of angles characterizing a given 
cone is four, precisely the number of degrees of freedom of 
the 2pcf. 


The 2pcf represented by this new set of angles can be 
expanded in the following way 

c(x,i, s, $) = e s (2i+i r +i V "' 

l 

x Pa(cosx)T>MM'(h,©,$)> (9) 

where D^ IM , is the Wigner D-matrix, Pi are the Leg¬ 
endre polynomials, and are the multipolar co¬ 

efficients of the expansion, which we term the angular- 
conic spectrum. We stress that the order of the angles 
( 77 , 0 , <f>) in D^ im , is important, since they are associ¬ 
ated to the eigenvalues M, L and M', in this order. 
Note also that the exchange symmetry Q now becomes 
C(x, V, ©: ( h) = C(Xi 77 ± 7 r, 0 , <h), which implies that M 
is even in the decomposition (|9|. 

In order for the multipolar coefficients C^ MM to be 
useful, they need to be related to the ai m defined in 
([ 2 ]), since the latter are more easily extracted from CMB 
maps. In principle, the relation among them can be found 
by equating Eq. © with Eq. © and by using the orthog¬ 
onality of the functions D^ IM , (rj, 0, <h) and P/(cosx) to 
write C ; LMM as a function of (ai irni ai 2m2 ). If performed 
naively, however, this task will lead to very complicated 
integrals coupling the conic angles (x,?7, 0,^) and the 
spherical angles (xi><AiiX 2 >^ 2 )- The easiest way to pro¬ 
ceed is to make use of the fact that the 2pcf is a scalar, 
and therefore the equality between Eqs. (J9| and ([2]) must 
hold in any coordinate system. We thus specialize the 
decompositions ([ 9 ]) and ([ 2 ]) to a coordinate system in 
which 0, $ and 77 are zero but x is not. Once the inte¬ 
gral over x is done, we rotate the system back to a general 
frame. This rotation can be performed using ( 77 , 0, $) as 
the three Euler angles, which can then be moved to the 
right-hand side of the equality using the orthogonality of 
the Wigner P-matrices. The details of this computation 
are given in the appendix. 
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The final result is 


-iLMM 


' =2n (“hr 




1 1 ,TO i 

h,'rri2 


L 

M' 


h 

mi 


h 

m 2 


Ji 


ILM 
I 1 I 2 ■ 


( 10 ) 

where Jj^ is a (non-square) matrix that couples the 
multipoles of the conic and spherical decompositions. 
This matrix comes entirely from the geometry of the 
problem, which means that its entries need to be com¬ 
puted only once. They are defined as 




ILM 

lih 


= £ 


L 

M 


h h 


V 

J h 


ml2Tn' 1 


(ii) 


P l {x)Pf i m {y)P% (: y)dx , 


where 

^l\ml2m' 'nii,—mVl2'rr, 

( 12 ) 

and y = \{x + l)/2] 1,/2 . The constant coefficients iy m are 
defined in ( |A1| ) . 

From equation ( flO] ) one sees that the expected values 
of the C^ MM relate linearly to the covariance matrix. 
Furthermore, in the limit of statistical isotropy, this re¬ 
lation becomes 


C, 


LMM' 


= CiSloSmoSm'o ■ 


(13) 


Similar to what happens with the BipoSH spectrum, 
any non-zero detection of the angular-conic spectrum 
with L > 0 is a measure of statistical anisotropy. How¬ 
ever, the multipole L has a simpler geometrical meaning, 
since it is associated with the orientation of the cone as 
defined in Fig. Q. The multipole L can be interpreted 
as an indication of a conic modulation of the 2 pcf over 
the CMB sky. This is an important aspect of the de¬ 
composition - that each of its angles has a clear geo¬ 
metrical meaning. One advantage is that it can make it 
simple to integrate one (or more) of them in a symmetric 
situation. The authors of Ref. [45], for example, con¬ 
sider a power-multipole test on the one-point function: 
C(n,h) = (AT 2 (h)). From the perspective of this work, 
their test is equivalent to Eq. (J9j) with X = y = 0. 

Finally, let us mention that the angular-conic spectrum 
is also linearly related to the BipoSH spectrum. Using 
Eqs. | 6 ]) and (10 1 , together with the transformation be¬ 
tween Clebsch-Gordan and Wigner 3J symbols (A4), we 
arrive at 


qLMM‘ 


' =2ttY 


(- 1 ) 


h-h-M 


I 1 J 2 


V2L + 1 




ILM 
1/2 ‘ 


(14) 


We will show in Sec. |III| how the angular-conic spectrum 
can be applied in a simple null test of statistical isotropy. 


C. Anisotropies through planar modulations 

Besides defining cones, two vectors n 1 and n 2 also de¬ 
fine a disc. The possibility of using a disc to represent 



Figure 2: Two unit vectors in the CMB sphere can be used to 
define a disc/plane. The number of angles characterizing this 
disc is three. Adding the isotropic angle \ between hi and 
h 2 , we arrive at the number of angles needed to characterize 
the most general 2pcf. 


the 2pcf was partially explored in Refs. |461I47| . Here we 
shall generalize these results to arrive at the most general 
2 pcf with planar symmetries. 

The geometry of the disc is characterized by three an¬ 
gles: two angles 0 and $ defining the overall orientation 
of the disc (i.e., its normal A 2 oc hi x n 2 ) and a third 
angle 77 measuring the rotation of the disc around its nor¬ 
mal 3 - see Fig. (pi. Including finally the angle \ between 
hi and n 2 , we have 

C (hr, h. 2 ) = C (x, 77 , 0 , 4>) , (15) 


where, again 

O<0<7T, 0<$<27T, 0<X< 7r > 0<77<27T. 


As previously done for the function ([7]) , the 2pcf func¬ 
tion defined by (15) can be expanded in terms of Wigner 
14-matrices and Legendre polynomials: 


C(x,?7,0,$) 


y- ( 2 £ + 1 )( 2 Z + l) LMM > 

l,m,m’ 1 

x P l {cosx)D% IM ,(r),e,$), (16) 


where are the multipolar coefficients of the ex¬ 

pansion. The exchange symmetry Q now becomes 
C(x, hi ©i $) = C(x,y,0 ± 77 , 4>), which further implies 
that in the decomposition (16 1 L must be even. The de¬ 
composition (16 1 generalizes the 2pcf introduced in Refs. 
|461 HS|, where the angle 77 was not included. 

Although the description so far seems to be identical 
with the conic representation of the 2 pcf, the relation 


3 We take r/ to be the angle from the great circle connecting the 
vectors z, —z and JV12, to the vector n\ along the disk. 
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between the disc’s angles (x, rj , 0 , $) and the spherical 
angles (xi> ^li X 2 i ^ 2 ) are different. Their interdepen¬ 
dence becomes clearer when expressing the new multipo¬ 
lar coefficients D ; iMM in terms of the covariance matrix 
(ai 

To find this relation we choose an initial coordi¬ 
nate system in which (x, 77 , 0 ,$) = (x, 0 , 0 , 0 ) and 
(Xi)^i)X 2 ,^ 2 ) = (?r/2, 0, 7 r/ 2 , fo). In this coordinate 
system the normal to the disc points in the z-direction 
and the remaining angle x is equal to </> 2 (see Eq. HD- 
After performing the integral over x, we rotate the sys¬ 
tem back to a general frame using ( 77 , 0 , <E>) as the three 
Euler angles, in that order (see the appendix for more 
details). This calculation gives 


M' 


__ ^l 2 m 2 ) 


L U 



where is a non-square matrix coupling the 

of angles involved. It is defined by 


tILM 

1 hh > 


(17) 
two set 


tILM _ ( L li l 2 \ A 

1 hh — / J \ M m m! ) ll iml 2 m' J 

m,m' ' ' 


(18) 


where 


h 1 ml 2 m’ — ^Zl m^l 2 


Pi{x)t 


im arccos x 


dx , 


i -1 


(19) 


In conclusion, the matrices jl L , AI and ll L , M can be seen 
as the weights that should be added to each eigenvalue 
of the BipoSH spectrum in order to obtain the angular- 
conic and angular-planar spectra, respectively. 


In what follows, it will be useful to introduce two new 
variables. Given the formal similarity between Eqs. (10 1 
and (Tt7|), we will define 


gLMM’ 


^ ] (^imd/ 2 ^ 2 ) 
l\,mi 
h,rri2 


L l\ I 2 
M’ mi m 2 


n lLM 

hh 


( 22 ) 

where Sf" AIAI stands for both the angular-conic and 
angular-planar power spectra, and Q^ AIM represents 
their respective coupling matrices. That is 


( rLMM' 
cLMM' _ I u l 


and 


n ILM _ 
VZiZ 2 - 


tILM 

J l\l2 
tILM 
, 1 hh 


(23) 


This will allow us to put all expressions in a unified 
description. 


Before we move on, let us illustrate the use of the 
angular-planar and angular-conic spectra to distinguish 
between different sources of anisotropy. First, let us con¬ 
sider anisotropic maps satisfying 

( a £i mi a C 2 m 2 ) = f{£\)f>hl 2 ±pb 771 1 777.2 (24) 


and A; m is a set of constant coefficients defined in (|A3|). 
Equation (fl7 


represents the desired relation between 
the multipolar coefficients in Eq. (16 ) and the temper¬ 
ature multipolar coefficients. In the limit of statistical 
isotropy these coefficients become 


where / is some predicted function of and p is any 
odd integer. This is a simple model of parity-violating 
anisotropy EL and can arise in many different theoret¬ 
ical contexts (T9i | 21 ] j49j [50]. For this covariance matrix 
the coefficients ([ 22 ]) become 


n^LMM' jr e e 

U x = L'IOlqOmoOM'O > 


( 20 ) 


as expected, since the multipole l measures the isotropic 
angular power of the CMB. 

We see here that, as it happens with the BipoSH and 
angular-conic spectra, the multipolar coefficients J)^ MM 
form a legitimate null estimator of statistical isotropy, 
since a measurement of any non-zero T)LAiAi ^ > 0 

is an indication of anisotropy. Given that l measures the 
isotropic angular power while L measures planar mod¬ 
ulations over an isotropic sky, the coefficients 

are called the angular-planar spectrum [4(7. Note that, 
again, the geometrical meaning of each angle involved in 
the construction of the 2 pcf is clear, allowing one to eas¬ 
ily marginalize over any desired degree of freedom in a 
symmetric situation. 

The angular-planar spectrum can be directly related to 
the BipoSH spectrum. Following the same computation 
leading to (141 we find 


V ^ MM ' =2ttJ2 — 


h-h-M 


li,h 


V2L+1 


-A, 


L-M 1 t ILM 


Ulo 


L hh 


( 21 ) 


(-'LMO 


5 


ILM 

v 


X 


£(-D 

_ 771 1 


l e 1 hTp\ 

0 ni\ ~mi J 


(25) 


Due to momentum conservation, the quantity inside 
square brackets is non-zero only when L is an odd in¬ 
teger. However, since the symmetry of the 2pcf restricts 
the planar multipole L to even values, for this particular 
example we have 

C ; LM0 ^ 0 and £>f M0 = 0 . (26) 


Likewise, models predicting a covariance matrix of the 
form {a eimi a} 2m2 ) = f(£ u such as hap¬ 

pens with CMB in the presence of a homogeneous mag¬ 
netic field EEL will lead to 

c lm, Tp = Q and v lm, Tp ^ Q _ ( 27 ) 

where we used the condition M' ± p = 0 imposed by the 
3J symbol. Since M' has to be an even number for the 
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angular-planar spectrum we conclude that Cf A/,=FP = 0 in 
this example. Evidently, exact results as above will not 
hold in practice, where all sorts of statistical noise and 
foregrounds might contribute differently to each spectra. 
For this one has to construct statistical estimators from 
the theoretical spectra which can be directly applied to 
a given CMB map. We next discuss how such estimators 
can be constructed. 


III. NULL TESTS OF ISOTROPY 

An interesting feature of the angular-conic and 
angular-planar spectra is that they can be used as null 
tests of isotropy that can potentially reveal the mecha¬ 
nisms producing the deviations from isotropy, thus giving 
hints on the mechanisms behind the observed signal. The 
geometrical interpretation of Eqs. ([7]) and (151 allows for 
the reduction of the number of angles in the 2pcf when¬ 
ever the peculiarities of the analysis permit. Most im¬ 
portant, though, is the fact that this feature allows us 
to control cosmic variance in a systematic way. Based 
on that, and having the reduction of cosmic variance in 
mind, in this work we will make the simplifying assump¬ 
tion that the angle ?? in Eqs. ([7| and (15) will not lead 
to significant modulations. In other words, we work with 


C(x, e, $) = E C(x, v, ©, $)dv , (28) 

which corresponds to taking the ry-monopole of Eqs. © 
and (16). Thus, from now on we shall use 


M = 0. 


(29) 


For convenience we can also drop the prime on M ', so we 
replace M' -¥ M, S) L0M -► Sf M and Q™ -»■ Q\^. 

Our primary motivation to assume Eq. (29) is sim¬ 


plicity, since it allows us to implement our method more 
easily. Nonetheless, it is important to justify this choice 
geometrically. Recall that, in the case of the planar 2pcf, 
the angle r) measures the rotation of the disc around its 
(fixed) normal. Thus, by assuming that M = 0 we will 
not be able to detect correlations of temperature along 
great circles in the sky, if they exist. In the case of the 
conic 2pcf the same angle will measure correlations of 
temperatures over small circular rings in the sky; again, 
such rings will not be detect if M = 0. Since it is not 
obvious that these features lie among known anomalies, 
this simplification seems appropriate in a first analysis. 
However, a thorough assessment of the CMB maps with 
the complete tools presented here can potentially reveal 
correlations of the type we are neglecting; the results 
of these analyses are in progress and shall be presented 
soon. 

A chi-square (null) test of conic/planar anisotropies 
can now be constructed. A simple unbiased estimator of 


the multipolar coefficients is 


gLM 


= E 


Q'l2lTl2 


li,mi 

^2,^2 


L l\ I 2 
M mi m 2 


Q\Th- (30) 


Its covariance around some expected theoretical value, 
is given by 

M LVMM‘ = ^ (§LM _ S L M y fgL'M' _ gL'M^ ^ 

Clearly, the most interesting theoretical model to test is 
ACDM, for which 


S, iM = 0, 


(L > 0) 


as follows from Eqs. (13) and (201. For this particu¬ 
lar model, and assuming Gaussianity of the temperature 
fluctuations, Eq. ([ 3 ]) holds. Then, with the help of Wick’s 
theorem, 


M 


ll' MM' 

l 


where 




8?r- 


2L + 1 


(of') Sll'Smm’, 

(31) 

E^^ (QwJ 2 - 

li,h 

(32) 


The fact that the matrix (311 is diagonal in the ACDM 


model is a consequence of the statistical independence of 
the ai m s in this model. As expected, this matrix depends 
exclusively on the angular power spectrum Cj, since this 
quantity completely defines the statistics in ACDM. 


Given the estimator (30) and its variance (311, we de¬ 
fine 


(x 2 )* - 


E 


gLM 


2L + 1 ^ ((jL\ 

M——L \ a l > 


(33) 


which is just a chi-square test divided by the 2L + l 
conic/planar degrees of freedom. Since by construction 

((x 2 )f) = 1, we define for simplicity 

(X 2 )N(X 2 )"-1- (34) 

Thus, any detection of (y 2 )^ yf 0 for L > 0 is an indica¬ 
tion of statistical anisotropy. 

An important remark is in order. If all the data one has 


is a single CMB map, the statistics (331 should be com¬ 


puted entirely in terms of that map’s data. In fact, this is 
the essence of a null test. Given a map, we treat it as if it 
were a ACDM map, and compute (33) accordingly. The 
computation of erf (or any other piece entering Eq. (33)) 
using a set of theoretical C)s, which supposedly generates 
the map at hand, will cease to be a null test, and will only 
bias our final result towards a priori expectations. Thus, 
in the case of a single map we compute with the power 
spectrum estimated by Ci = (21 + l) -1 )T) m \ai m \ 2 . 
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Figure 3: Angular-conic null test of isotropy applied to the 
Planck Commander 2015 temperature map, along with its 
mask. The plot shows the quantity (\ 2 )f versus L for dif¬ 
ferent angular multipoles l. The shaded contours represent 
2 a cosmic variance from 10 3 FFP6 simulations using the same 
mask. We also performed analyses with the other three Planck 
foreground-cleaned 2015 maps: SMICA, NILC, and SEVEM, ob¬ 
taining qualitatively similar results. 


A. Null tests of isotropy with Planck data 

In this section we perform a simple statistical analy¬ 
sis of the 2015 Planck data release [2] using the tools we 
developed. While this analysis is not intended to be ex¬ 
haustive, it offers an important sanity check of the avail¬ 
able data. For these analyses we have used the Planck 
Commander 2015 temperature map, along with its mask. 
The data points were compared with the 1000 Full Focal 
Plane (FFP6) simulations provided by the Planck team, 
to which we have also applied the Commander-Ruler 
mask, so as to ensure that we are comparing quantities 
with the same foreground treatments. The results of our 
analysis are shown in Figs. © and Q, where the data 
points are compared to the 2cr cosmic variance bars from 
the FFP6 simulations. We have also performed analy¬ 
ses with the other three Planck foreground-cleaned 2015 
maps: SMICA, NILC, and SEVEM, obtaining qualitatively 
similar results. 

Our results show no drastic discrepancies between the 
Planck 2015 data and the ACDM model at the multi¬ 
poles we tested, although it is interesting to note that 
the cases {x 2 )iL 3 for the cone and ( X 2 )b =4 5 f° r the disc 
present a consistent deficit of conic/planar modulations 
in the scales we considered. It is important to mention 
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L L 

Figure 4: Angular-planar null test of isotropy applied to 
the Planck Commander 2015 temperature map, along with its 
mask. The plot shows the quantity (y 2 )f versus L for differ¬ 
ent angular multipoles l. The shaded contours represent 2a 
cosmic variance from 10 3 FFP6 simulations using the same 
mask. We also performed analyses with the other three Planck 
foreground-cleaned 2015 maps: SMICA, NILC, and SEVEM, ob¬ 
taining qualitatively similar results. Note that, for this test, 
the planar multipole L is always even. 

that, while the conic/planar multipoles L are statisti¬ 
cally independent in the ACDM model (see Eq. ©), 
the angular multipoles l are not. In other words, while 
the points in each frame of Figs. ([3| and Q are inde¬ 
pendent from each other, they are not independent from 
the points in other frames of the same figure. Given that 
there are 16 independent points in the multipole range 
we considered, on average, only 16 x 4.55% = 0.73 of 
them should fall outside the 2a (95.45% Cl) variance er¬ 
ror bars. This analysis is in clear agreement with our 
findings. 


IV. CONCLUSIONS AND PERSPECTIVES 

The impressive agreement of the ACDM model with 
current CMB data compel us to go beyond the simple 
statistical framework of a Gaussian, homogeneous and 
isotropic universe. Assuming that CMB is Gaussian, a 
fact which is supported by the latest Planck data, hints 
of new physics can be found in the realm of statistical 
anisotropy - a possibility which is still open to debate. 

In this work we have proposed two new representa¬ 
tions of the 2pcf as alternatives to the popular bipolar 
power spectrum (BipoSH) analysis. Besides being model 






































independent, these tools are entirely based on the geom¬ 
etry defined by the two unit vectors in the CMB sphere, 
namely, a cone and a disc. These tools differ from the 
BipoSH analysis in two main aspects. First, the new 
decompositions of the 2 pcf are geometrically inspired, 
which means that null tests of isotropy based on their 
multipolar coefficients can help to elucidate the physical 
mechanism behind signals of anisotropy. We have illus¬ 
trated this feature with concrete examples of anisotropic 
maps in which only one of these estimators will be non¬ 
zero. Thus, for example, a statistically significant detec¬ 
tion of the angular-planar spectrum cannot result from 
parity-violating physics, since in this case the angular- 
planar spectrum is zero for all multipoles. Second, the 
clear geometrical role of the anisotropic degrees of free¬ 
dom used as variables in the 2 pcf allows us to construct 
simpler correlation functions whenever we have a more 
symmetric situation at hand. This feature has the impor¬ 
tant consequence of allowing us to reduce cosmic variance 
in a systematic way. 

The angular-conic and angular-planar null tests of 
isotropy have not shown significant deviations of sta¬ 
tistical isotropy in the lowest multipoles of the Planck 
data, although some angular multipoles presented inter¬ 
esting low conic/planar modulations. Our analysis does 
not reveal any new statistical anomaly, although known 
features, such as the hemispherical power anomaly, are 
expected to appear at higher angular multipoles , which 
were not included in our first analysis. 

A more complete analysis of existing observational 
data, including a larger range of multipoles and a bet¬ 
ter assessment of systematics, is postponed to a future 
work. Of particular interest is the application of the 
angular-conic spectra in the investigation of the hemi¬ 
spherical power asymmetry found in Refs. (52W54j . In¬ 
deed, these references considered a pixel-based variance 
estimator which resembles in many ways the conic de¬ 
grees of freedom that we have introduced. In this re¬ 
spect, it is worth mentioning that the conic 2 pcf might 
have some relevance for the detection of Baryonic Acous¬ 
tic Oscillations. Indeed, the ring-like pattern that BAO 
produces in the distribution of galaxies leads exactly to 
a three-dimensional cone centred on the observer. 
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Appendix A 


We present here the derivation of our main results, 
Eqs. (10) and (Jl7|. We also collect useful identities and 
mathematical formulas used in the text. 


1. Spherical Harmonics and Wigner 3J symbols 


Our definition for the spherical harmonics is 

Y lm (xA) = {-l) m rnmPr^os X )e im \ 

where 

121 + 1 (l — to)! (I — to)! 


V Im — 


47 t (l + m)\ m (l + to)! 

At the point (y, <f>) = ( 7 r/2, </>), it simplifies to 
Y lm (TT/ 2 ,ct>) = x lm e im<l> , 

where 


(Al) 


(A2) 


, s. ( f + m ) 

A lm = (“I) 2 X 


2 i+l (I— 

4 -k (Z+m)!! (/ —m)!! 


1/2 


l + m = even 
otherwise. 


(A3) 


The relation between Wigner 3Js and Clebsch-Gordan 
coefficients are 


h h h 

m i TO 2 — m 3 


(- 1 ) 


h—h—ms 


_ £yl3 rn 3 


sJ 21 3 + 1 hm.hm. 

Other useful identities include 

^ = V2iTTs L0 , 

(-i y~ m 




L l l 

0 TO —TO 


0 l l 

0 TO TO J y/21 + l ' 


(A4) 


(AS) 


(A 6 ) 


We also remind two useful orthogonality relations of 
the Wigner rotation matrices. These are 


and 


f D l1 , ( w)D l2 , (w)D l3 , (w)dw = 

I mimj A / 1712^2 ' m^m 3 V / 


87T 


2 I h h h \ ( h h h 


m'i m' 2 TO 3 ) \ mi m 2 m 3 


with w = (a, p, 7 ) being the three Euler angles. 















9 


2. Derivation of Eq. (10) for conic-angular 
modulations 


After expanding the 2pcf in terms of Legendre polyno¬ 
mials and Wigner D-matrices, we equate expressions Q 
and ([2]). Since the 2pcf is a scalar, this equality should 
hold in a coordinate system in which the symmetry axis 
of the cone is aligned with the 3 -axis. That will mean: 


Combining this result with (fA 8 |, we finally get (13). 

3. Derivation of Eq. ( |17| ) for planar-angular 
modulations 


Equating the expanded 2pcf ( fl 6 | with ([2]), we choose a 
particular coordinate system in which the normal to the 
disc points in the 3-direction, which means: 


0 = $ = 77 = 0 , xi = X2 = X / 2 , fa = fa ~ 7T = 0 . 

Using the identity D^ M , (0,0,0) = Smm' and the or¬ 
thogonality of the Legendre polynomials, we then arrive 
at 


0 = $ = n = 0, Xi = X 2 = tt/2, (pi =0, X = fa¬ 
ll sing 0 , 0 , 0 ) = S M 'm, Y lm fa/2,(p) = A 

and integrating over P;(cosy), we find 


E 

L,M 


(2L + 1) 

27r 


qLMM 


'y ' ( a l 1 m 1 a l 2 m 2 ) Jl 1 m 1 l 2 m 2 1 

h ,m 1 
h,fri2 


_ (A7) 

where j] m m was defined in Eq. (12). If we now 
rotate the axes back to a general coordinate system using 
u> = { 77 , 0, $} as the three Euler angles, the coefficients 
qLMM anc j aim w j]j transform as [38] 


E 

L,M 


(2L + l) 

27r 


sj-yLM M 


E 17711^/2^2) Ilimil 2 rri 2 

li ,?ni 

Z2,Wl2 


where l\ m i m was introduced in (19). We now rotate 
back to a general coordinate system using w = { 77 , 0 , $} 
as the Euler angles and the fact that, under rotations 




LMM 1 t~\L ( \ 




>=E 


O-lm'Dram 1 ( w ) • 


M' 


Then, we multiply both sides of (A7) by and 

use the orthogonality of the Wigner rotation matrices to 
isolate C^ MM . Finally, we use (—1) M+M qL,-m-m _ 
qLMM *, w } 1 i c ] 1 follows from the reality of the 2 pcf, and 
substitute —> C^ MM ', since we could equally well 

have started with the complex conjugate in the expansion 
©■ Dropping primes and tildes, we arrive at (|10|). 


M' 771' 

From this point on, the deduction is similar to the case 
of the cone. After some redefinitions and relabeling of 
the indices, we finally arrive at (|l7|). 


a. Isotropic limit 


a. Isotropic limit 

In order to derive the isotropic limit ( |l3| ), we note that, 
in this limit, (u^ imi <z^ 2m2 ) — ( 1 ) 2 C 
Since the Wigner 3J symbol appearing in Eq. (10) is zero 


unless M' + mi + m 2 = 0, this implies that M' = 0. 
Then, using Eq. (]A5[), we arrive at 


c lmm' = 2irJ2Ci 1 (-l) h V2h+ij!^ x S LO S M 'o ■ 


In the above expression, the non-vanishing terms of 
the coupling matrix are of the form ■ Combining 
Eq. (A 6 ) with the addition theorem for the associated 


Legendre polynomials one can show that 

1 


J “" = Pmfafa s,hSM °- 


The derivation of (|20j) follows a similar deduc¬ 
tion to the one of the angular-conic power spec¬ 
trum. In this limit, the expression (ae irni ae 2 m 2 ) = 
(—1 ) m2 Ce 1 6e 1 t 2 6 rni - rn2 implies that M' = 0. Then, us¬ 
ing Eq. (A51, we arrive at 


p lmm' = 27 r^C' il (-iyV2/i + l/^f x 6 lo S M 'o ■ 

(A10) 


(A 8 ) Using again Eq. (A 6 |, it is not difficult to show that 


tIOM _ 
1 hh - 


(- 1)* 1 1 
y/2li + 1 27T 


Slh^MO ■ 


(All) 


(A9) Combining this result with (AlOl, we finally get (20 1 . 
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